Multifractal Structure of the Harmonic Measure of Diffusion Limited Aggregates 
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The method of iterated conformal maps allows to study the harmonic measure of Diffusion Lim- 
ited Aggregates with unprecedented accuracy. We employ this method to explore the multifractal 
properties of the measure, including the scaling of the measure in the deepest fjords that were hith- 
erto screened away from any numerical probing. We resolve probabilities as small as 10 -35 , and 
present an accurate determination of the generalized dimensions and the spectrum of singularities. 
We show that the generalized dimensions D q are infinite for q < q* , where q* is of the order of 
—0.2. In the language of /(a) this means that a max is finite. The f(a) curve loses analyticity 
(the phenomenon of "phase transition") at a m ax and a finite value of f(a maK ). We consider the 
geometric structure of the regions that support the lowest parts of the harmonic measure, and thus 
offer an explanation for the phase transition, rationalizing the value of q* and /(a max ). We thus 
offer a satisfactory physical picture of the scaling properties of this multifractal measure. 



I. INTRODUCTION 

Multifractal measures are normalized distributions ly- 
ing upon fractal sets. Such measures appear naturally in 
a variety of nonlinear-physics context, the most well stud- 
ied being natural measures of chaotic dynamcial systems 
|], ||, ||. Other well studied examples are the voltage 
distribution of random resistor networks |5| . In this 
paper we address the harmonic measure of Diffusion Lim- 
ited Aggregates Q , which is the probability measure for 
a random walker coming from infinity to hit the bound- 
ary of the fractal cluster. This was one of the earliest 
multifractal measures to be studied in the physics liter- 
ature , but the elucidation of its properties was made 
difficult by the extreme variation of the probability to 
hit the tips of a DLA versus hitting the deep fjords. 
With usual numerical techniques it is quite impossible 
to estimate accurately the extremely small probabilities 
to penetrate the fjords. Contrary to harmonic measures 
of conformally invariant fractals like random walks and 
percolation clusters whose multifractal properties can be 
solved exactly ||, E|, the present multifractal measure 
posed stubborn barriers to mathematical progress. 

The multifractal properties of fractal measures in gen- 
eral, and of the harmonic measure of DLA in particular, 
are conveniently studied in the context of the generalized 
dimensions D q , and the associated /(a) function |n], [Ll| . 
The simplest definition of the generalized dimensions is 
in terms of a uniform covering of the boundary of a DLA 
cluster with boxes of size i, and measuring the probabil- 
ity for a random walker coming from infinity to hit a piece 
of boundary which belongs to the i'th box. Denoting this 
probability by Pi(£), one considers |l(| 



D q = lim 



i logEi^W 



e->o q - 1 



loge 



(i) 



where the index i runs over all the boxes that contain 
a piece of the boundary. The limit Dq = lim q _ t0 + D q 
is the fractal, or box dimension of the cluster. D\ = 
linig— >i+ D„ and D2 are the well known information and 



correlation dimensions respectively ^ [l2|, [ll|. It is well 
established by now jllj that the existence of an interest- 
ing spectrum of values D q is related to the probabilities 
Pi(£) having a spectrum of "singularities" in the sense 
that Pi(£) ~ £ a with a taking on values from a range 
amin < ol < a max . The frequency of observation of a 
particular value of a is determined by the the function 
f(a) where (with r(q) = (q — l)D q ) 



f{a) =aq(a) -T[q(a) 



dr(q) 
dq 



a(q) . (2) 



The understanding of the multifractal properties and 
the associated f(a) spectrum of DLA clusters have been 
a long standing issue. Of particular interest are the val- 
ues of the minimal and maximal values, a m i n and a maxi 
relating to the largest and smallest growth probabilities, 
respectively. 

The minimal value of a is relatively easy to estimate, 
since it has to do with the scaling of the harmonic mea- 
sure near the most probable tip. While the often cited 
Turkevich-Scher conjecture ]l4| that a m m satisfies the 
scaling relation Do = 1 + a m - ln is probably not exact, it 
comes rather close to the mark. On the other hand, the 
maximal value of a is a much more subtle issue. As a 
DLA cluster grows the large branches screen the deep 
fjords more and more and the probability for a random 
walker to get into these fjords (say around the seed of the 
cluster) becomes smaller and smaller. A small growth 
probability corresponds to a large value of a. Previ- 
ous literature hardly agrees about the actual value of 
&max- Ensemble averages of the harmonic measure of 
DLA clusters indicated a rather large value of a max ~ 8 
|l5| . In subsequent experiments on non-Newtonian flu- 



and on viscous fingers 1 7 , similar large values 



ids ]1 

of Oi max were also observed. These numerical and exper- 
imental indications of a very large value of (x max led to a 
conjecture that, in the limit of a large, self-similar clus- 
ter some fjords will be exponentially screened and thus 
causing a max — ► 00 |Jj . 

If indeed a max — > 00, this can be interpreted as a phase 
transition [fij| (non-analyticity) in the q dependence of 
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D q , at a value of q* satisfying q* > 0. If the transition 
takes place for a value q* < then a max is finite. Lee and 
Stanley (20) proposed that a max diverges like R 2 /\tlR 
with R being the radius of the cluster. Schwarzer et al. 
pl| proposed that a max diverges only logarithmically in 
the number of added particles. Blumenfcld and Aharony 
p2[ proposed that channel-shaped fjords are important 
and proposed that a max ~ in m wnere M is the mass 
of the cluster; Harris and Cohen p3[ , on the other hand, 
argued that straight channels might be so rare that they 
do not make a noticeable contribution, and a max is finite, 
in agreement with Ball and Blumenfeld who proposed 
p4j that a max is bounded. Obviously, the issue was not 
quite settled. The difficulty is that it is very hard to 
estimate the smallest growth probabilities using models 
or direct numerical simulations. 

In a recent Letter we used the method of iterated 
conformal maps to offer an accurate determination of the 
probability for the rarest events. The main result that 
was announced was that a max exists and the phase tran- 
sition occurs at a q value that is slightly negative. In the 
present paper we discuss the results in greater detail, and 
offer additional insights to the geometric interpretation 
of the phase transition. In Sect. 2 we summarize briefly 
the method of iterated conformal maps and explain how 
it is employed to compute the harmonic measure of DLA 
with unprecedented accuracy. In Sect. 3 we perform the 
multifractal analysis and present the calculation of the 
f(a) curve. In Sect. 4 we discuss a complementary point 
of view of the scaling properties of the rarest regions of 
the measure, to achieve in Sect. 5 a geometric interpre- 
tation of the phase transition. Sect. 6 offers a short 
discussion. 



done, but nevertheless this is the basic principle of our 
approach. 

We thus need to find the conformal map 
excellent method for this purpose was developed in a re- 
cent series of papers |^6|, |2?], ^8). The map $( n )(tt>) is 
made from compositions of elementary maps 4>\fii 

$(»>(«,) =$(»-D (0 a „, 9t ») , (4) 

where the elementary map cf>x.e transforms the unit circle 
to a circle with a semi-circular "bump" of linear size Vx 
around the point w — e % ® . We use below the same map 
d>\ a that was employed in @ || [29[ |(} . With this 
map $(")(io) adds on a semi-circular new bump to the 
image of the unit circle under $( n-1 ) (w). The bumps in 
the z-plane simulate the accreted particles in the physical 
space formulation of the growth process. Since we want 
to have fixed size bumps in the physical space, say of 
fixed area Ao, we choose in the nth step 

A = A ° (51 

n |$(»-i)'(e i9 ")| 2 ' 

The recursive dynamics can be represented as iterations 
of the map <f>\ n ,e n (w), 

$ (n) (w) = Al ,0i o A2:e2 o . . . o (j) Xnfi Juj) . (6) 

It had been demonstrated before that this method repre- 
sents DLA accurately, providing many analytic insights 
that are not available otherwise j2{| |3(| . 

B. Computing the harmonic measure 



II. ACCURATE CALCULATION OF THE 
HARMONIC MEASURE 

A. DLA via iterated conformal maps 

Consider a DLA of n particles and denote the bound- 
ary of the cluster by z(s) where s is an arc-length param- 
terization. Invoke now a conformal map which 
maps the exterior of the unit circle in the mathematical 
plane ui onto the complement of the cluser of n particle in 
the z plane. On the unit circle e %e the harmonic measure 
is uniform, P(8)d0 = d9/2ir. The harmonic measure of 
an element ds on the cluster in the physical space is then 
determined as 



P(s)ds 



ds 



|$'(™)(e ie )| 



(3) 



where $( n )(e ze ) = z(s). Note that in electrostatic par- 
lance 1/1$ ( n )(w)| is the electric field at the position 
z = $( Ti )(a;). Thus in principle, if we can have an ac- 
curate value of the conformal map (w) for all values 
ui — e l6 we can compute the harmonic measure with de- 
sired precision. We will see that this is easier said than 



In terms of computing the harmonic measure we note 
the close relationship between Eqs.(||) and (||). Clearly, 
moments of the harmonic measure can be computed from 
moments of X n . For our purposes here we quote a result 
established in |27| , which is 



(XD = (1/2tt) X n (0)dd~n 
Jo 



-2qD 2q+1 /D 



To compute r(g) we rewrite this average as 

A« +1 /2( s ) 



(Xn 



1\ - 



ds 



XUs) = ds 



(7) 



(8) 



where s is the arc-length of the physical boundary of 
the cluster. In th e last equality we used the fact that 
\dd/ds\ — \J X n /Xo. We stress at this point that in order 
to measure these moments for q < we must go into 
arc-length representation. 

To make this crucial point clear we discuss briefly what 
happens if one attempts to compute the moments from 
the definition (Q). Having at hand the conformal map 
§( n )(e l6 ), one can choose randomly as many points on 
the unit circle [0, 2it] as one wishes, obtain as many (ac- 
curate) values of A„, and try to compute the integral 
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FIG. 1: panel a: the boundary of the cluster probed by a 
random search with respect to the harmonic measure. Panel 
b: the boundary of the cluster probed by the present method. 



as a finite sum. The problem is of course that using 
such an approach the fjords are not resolved. To see this 
we show in Fig.l panel a the region of a typical clus- 
ter of 50 000 particles that is being visited by a random 
search on the unit circle, with 50 000 samples. Like in 
direct simulations using random walks, the rarest events 
are not probed, and no serious conclusion regarding the 
phase transition is possible. Another method that can- 
not work is to try to compute (A^) by sampling on the 
arc-length in a naive way. The reason is that the inverse 
map -1 (s) cannot resolve 9 values that belong to 

deep fjords. As the growth proceeds, reparametrization 
squeezes the 9 values that map to fjords into minute in- 
tervals, below the computer numerical resolution. To see 
this, consider the following estimate of the resolution we 
can achieve in the physical space, 



AO = 



A& 



or equivalently 



$(«)' 



As 



= VXr. 



A9 



As 

Va 



(9) 



(10) 



FIG. 2: A typical growth process in which an existing branch- 
cut is "burried" under the new bump. Such events reduce the 
number of branch cuts below 2n with n being the number of 
particles. 



On the left hand side we have the resolution in the phys- 
ical space relative to the fixed linear size of the particles. 
With double precission numerics we can resolve values of 
A9 ~ 10 -16 and since we know that the values of A30000 
can be as small as 10~ 70 inside the deepest fjords (and 
see below), we see that 
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As 



Va 10- 



-16 



10 



1!) 



(11) 



Therefore the resolution in the physcial space which is 
necessary to achieve a meaningful probe of the deep fjord 
is highly inappropriate. 

The bottom line is that to compute the values of A n (s) 
effectively we must use the full power of our iterated con- 
formal dynamics, carrying the history with us, to iterate 
forward and backward at will to resolve accurately the 9 
and X n values associated with any given particle on the 
fully grown cluster. 

To do this we recognize that every time we grow a semi- 
circular bump we generate two new branch-cuts in the 
map . We find the position on the boundary between 
every two branch-cuts, and there compute the value of 
A„. The first step in our algorithm is to generate the loca- 
tion of these points intermediate to the branch-cuts |H] . 
Each branch-cut has a preimage on the unit circle which 
will be indexed with 3 indices, wA = e^p[i6A']. The 
index j represents the generation when the branch-cut 
was created (i.e. when the jth particle was grown). The 
index I stands for the generation at which the analysis is 
being done (i.e. when the cluster has £ particles). The in- 
dex k represents the position of the branch-cut along the 
arc-length, and it is a function of the generation I. Note 
that since bumps may overlap during growth, branch- 
cuts are then covered, cf. Fig. 0. Therefore the maximal 
k, k max < 11. After each iteration the preimage of each 
branch-cut moves on the unit circle, but its physical po- 
sition remains. This leads to the equation that relates 
the indices of a still exposed branch-cut that was created 
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at generation j to a later generation n: 

= $(«) o . . . o ^ li9 , +1 («^)) 

- . (12) 



Note that the sorting indices k(j) are not simply related 
to k(n), and need to be tracked as follows. Suppose 

that the list w^Zp is available. In the nth genera- 
tion we choose randomly a new 9 n , and find two new 
branch-cuts which on the unit circle are at angles 0^. If 
one (or very rarely more) branch-cut of the updated list 

(j)^ 1 g (uij^n-P) is covered, it is eliminated from the list, 
and together with the sorted new pair we make the list 

k(n) 
W ■ . 

J,n 

Next let us find the mid-positions at which we want to 
compute the value of A ra . Having a cluster of n particles 
we now consider all neighboring pairs of preimages w,- „ 



and ' 



k(n) + l 



, that very well may have been created at two 



different generations j and J. The larger of these indices 
( J without loss of generality) determines the generation 
of the intermediate position at which we want to compute 

the field. We want to find the preimage uf^ of this 
mid-point on the unit circle , to compute Xk( n ) there 
accurately. Using definition (fL2) we find the preimage 



/ k(n) 

ar §( w j,j 



[arg(^ ( / ) )+arg(^ ( J /)+1 )]/2. (13) 



In Fig.l panel b we show, for the same cluster of 50 000, 
the map <$>( J \ujj ) with k(n) running between 1 and 
fcmax, with J being the corresponding generation of cre- 
ation of the mid point. We see that now all the particles 
are probed, and every single value of Xu n ) can be com- 
puted. 

To compute these Xun) accurately, we define (in anal- 
ogy to Eq.(12)) for every J < m < n 



k(n) 
L J,rn 



(14) 



Finally Xk( n ) is computed from the definition (^j) with 

V \ U .J,n ) — 9\ n> 8 n \ U J,n ) ■ ■ ■<P\ J + 1 ,9j +1 \ u J,J+l) 



k(n) 
J,J 



(15) 



We wish to emphasize the relevance of this equation: 
the problem with the coarse resolution that was exposed 
by Eq. (11) occurs only inside the deepest fjords. We 
note however that the particles inside the deep fjords 
were deposited when the clusters were still very small. 
For small clusters the resolution of the fjords does not 
pose a difficult problem. Therefore, when we evaluate 
the derivative $^ n ^ inside the deepest fjord at a point 
we make use of the fact that J<n and write the 



u 



k(n) 
J.n 



derivative in the form 

$ (»)' ( u Hn)^ = $ (J)' £ u fc(n)j x " correcting tcrms » ( 16 ) 



On the left hand side of Eq. (pj) we see that within our 
limited numerical resolution u 



n) fc(n)+l 



J.n 



a 



J.n 



and the cor- 

rcponding values of A„ are almost identical whereas for 
the righthand side this is not the case. By keeping track 
of the branch cuts we improve the precision inside the 
fjords dramatically. In other words, the large screening 
inside the fjords is simultaneously the problem and the 
solution. The problem is that we cannot use the stan- 
dard approach in evaluating $W . The solution is that 
for a point inside the deepest fjords we always have 



that J <§; n and therefore the evaluation ([L5|) helps to 
improve the resolution. 

In summary, the calculation is optimally accurate since 
we avoid as much as possible the effects of the rapid 
shrinking of low probability regions on the unit circle. 
Each derivative in ( |l5| ) is computed using information 
from a generation in which points on the unit circle are 
optimally resolved. 

The integral (||) is then estimated as the finite sum 
V^o J2k(n) ^k(n) ■ W e should stress that for clusters of the 
order of 30 000 particles we already compute, using this 
algorithm, Afe( n ) values of the order of 10~ 70 . To find the 
equivalent small probabilities using random walks would 
require about 10 35 attempts to see them just once. This 
is of course impossible, explaining the lasting confusion 
about the issue of the phase transition in this problem. 
This also means that all the f(a) curves that were com- 
puted before |l5|, |32]] did not converge. Note that in our 



calculation the small values of A 



k{n) 



are obtained from 



multiplications rather than additions, and therefore can 
be trusted. 



III. MULTIFRACAL ANALYSIS OF THE 
HARMONINC MEASURE 

Having the accurate values Afc( n ) we can now compute 
the moments (Q). Since the scaling form on the RHS 
includes unknown coefficients, we compute the values of 
r(q) by dividing (A%) by (A|), estimating 



r(q) » -D 



log(A£)-log(A? 
log n — log n 



(17) 



Results for r(q) for increasing values of n and n are shown 
in Fig. ||, panel a. It is seen that the value of r(g) appears 
to grow without bound for q negative. The existence of 
a phase transition is however best indicated by measur- 
ing the derivatives of r(q) with respect to q. In Fig. ^ 
panel b we show the second derivative, indicating a phase 
transition at a value of q that recedes away from q = 
when n increases. Due to the high accuracy of our mea- 
surement of A we can estimate already with clusters as 
small as 20-30 000 the q value of the phase transition as 
q* = —0.18 ± 0.04. It is quite possible that larger clus- 
ters would have indicated slighly more negative values 
of q* , (and see below the results of different methods of 
estimates), but we believe that this value is close to con- 
vergence. The fact that this is so can be seen from the 
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FIG. 3: Panel a: the calculated function r(q) for clusters of 
n particles, with n =10 000,15 000, 25 000 and 30 000. Panel 
b: the second derivative of r(q) with respect to q. 



FIG. 4: The calculated function /(a) using r(q) calculated 
from a cluster with n — 30 000 particles. This f{a) is almost 
indistinguishable from the one computed with n = 25000 par- 
ticles. We propose that this function is well converged. The 
black dot denotes where the curve ends, being tangent to the 
line with slope -0.18. 



IV. ALTERNATIVE WAY TO APPROACH THE 
PHASE TRANSITION 

An alternative way to the multifractal analysis is ob- 
tained by first reordering all the computed values of \k( n ) 
in ascending order. In other words, we write them as a 
sequence {A„(i)}i e j where / is an ordering of the indices 
such that \ n (i) < X n (j) if i < j. The number of samples 
we consider is usually large and therefore the discrete in- 
dex i/N might be treated as a continuous index < x < 1 
and X n as a non-decreasing function of x, 

A„ = f{x) (18) 



f(a) curve which is plotted in Fig. |j. A test of conver- 
gence is that the slope of this function where it becomes 
essentially linear must agree with the q value of the phase 
transition. The straight line shown in Fig. 3 has the slope 
of -0.18, and it indeed approximates very accurately the 
slope of the f(a) curve where it stops being analytic. The 
reader should also note that the peak of the curve agrees 
with D ps 1.71, as well as the fact that r(3) is also D as 
expected in this problem. The value of a max is close to 
20, which is higher than anything predicted before. It is 
nevertheless finite. We believe that this function is well 
converged, in contradistinction with past calculations. 



We next consider the distribution of p(X n ), which is 
calculated by the usual transformation formula 

p(X n ) ~ J 5{\ n - f{x))dx = yjT^jyf 

Using the distribution of A„ we now do the following 
rewritings, 

/ \ q n d6 = / X q n — ds~ / Xf 2 ds (20) 
Jo Jo as J Q 




Xr*p(X n )dX n . (21) 



In Fig. g we will show that our function f(x) obeys a 
power law for low values of x, 

f(x) ~ x 13 , for x < 1 . (22) 
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FIG. 5: Values of An sorted in ascending order with respect 
to the variable x = i/N. This function is a pure power law 
for values of log x smaller than the position of the circle. The 
power law is characterized by an exponent f3 ~ 8.55. This is 
consistent with a phase transition for q* ~ —.23 



This in turn implies a power law dependence of p(X n ) on 

An 

p(\ n ) — ~ , for A n « 1 . (23) 

(z(An)r 

This power law tail means that the moment integral 
( pl| ) diverges for values of q below a critical value q c given 

by, 

+ 1 (24) 

Thus q c = — \ — 4. From (0) we see that the value of q* 
satisfied the relation 

q* = 2q c + 1 = -2/P . (25) 

In Fig. H we show how the values of A„ depend on x for 
small values of x. The data is taken from a cluster with 
n = 20000. Denoting the value of A„ that is marked as a 
full circle by A c , the figure supports the existence of the 
power law ( p2| ) for values of A„ smaller than A c . Needless 
to say this alos implies that p(X n ) scales according to 
Eq. (p3[) . By averaging over 16 clusters of size n=20000 
we estimate the slope in Fig. || to be m 8.55 or 

q* = -0.23 ± .05 . (26) 

Obviously, this result is in agreement with our direct cal- 
culation in Sect. 3. 



V. GEOMETRICAL INTERPRETATION OF 
THE PHASE TRANSITION 

At this point we would like to interpret the origin of 
the phase transition, which in light of the last section 
stems from the power law behavior of p(X n ) for small 
values of A„. We first identify the region on the cluster 
that supports the low values of A„ that belong to the 
power-law tail of p(A„). 

Consider again Fig. |[ The point denoted above as 
A c defines the maximum value for which we see a power 
law in A„ vs. x. Therefore the set responsible for the 
phase transition is the union of bumps with a value of 
A„ for which A„ < A c . This set is referred to below as 
the "critical set" , and is shown in figure q, both on the 
background of the rest of the cluster, and as an isolated 
set. This figure suggests a geometric interpretation: the 
fjords in the figure all seem to have a characteristic angle. 
We will try first to confirm this impression using careful 
numerics. 

Clearly, the set has several fjords; we consider them 
individually. Fig. |?j shows an example of such a fjord. 
For each fjord we find the point with the minimum prob- 
ability and use this for defining the bottom (or deepest 
point). Second, from the inside, we move to the two adja- 
cent points which together with the deepest point define 
an angle. This angle is recorded, and we move to the 
next pair of points, and so on until the value of A„ ex- 
ceeds A c . Fig.^| panel b shows how the angle varies with 
the number of steps k. For most of the fjords considered 
the angle is quite large for a small number of steps (up 
to 3-4 steps). As more steps are taken the angle settles 
on a characteristic value around which it fluctuates. For 
a larger number of steps we reach the outer parts of the 
fjord and the angle does no longer reflect the geometry 
inside the fjord. The dependence of the angle on k as 
shown in Fig. |?] is typical for all the fjords of the set 
causing the phase transition and therefore we see a peak 
in the distribution of all the measured angles. This peak 
identifies an angle which is characteristic to the fjords. 
Fig. U is the distribution of such typical angles over one 
cluster. We determine the characteristic angle, say j c by 
locating the maximum of the distribution. Finally, we 
calculated the average of the charateristic angle 7 C over 
15 clusters of size n — 20000. On the basis of that we 
determine the angle to be 

7c = 35° ± 6° . (27) 

Finally, we can offer a geometrical model to interpret 
the phase transition. The results presented in this sec- 
tion indicate that to a reasonable approximation the least 
accessible fjords can be modeled as a wedge of included 
angle 7 C . In Appendix [A] we compute the power law ex- 
pected for p(A„) for a wedge. The final result is 



p(X n ) ~ A„ 2(Q ~ 1) , (28) 
where 7 C = ir/a. Using our numerical result for j c and 
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FIG. 6: The set of all particles which are associated with 
values of A n belonging to the power law region shown in Fig. 
|E| In panel a we show the set on the background of the cluster 
in grey, and in panel b the set isolated from the rest of the 
cluster. 



Eq.@ we predict finally q* = -.24 ± .05. Obviously, 
this is in excellent agreement with the previous findings. 

In addition, we should comment on the interpretation 
of f{a max ) which is the value of the f(a) curve at the 
point of loss of analyticity. Within the wedge model of- 
fered here, this must be the fractal dimension of the set 
of wedges that support the scaling law (p8|). We have 
attempted to determine this dimension numerically by 
counting the number of fjords in the critical set shown in 
Fig. (0) as a function of the number of paricles n in the 
cluster. While the result of such a calculation is consis- 
tent with the proposition, the available statistics is not 
sufficient to establish it firmly. We thus conclude with 
the proposition as a conjecture, i.e. that f(a max ) can 
be interpreted as the dimension of the set of fjords that 



FIG. 7: Panel a: A typical deep fjord resolved on the sacle of 
the particles. From the deepest particle the angle is computed 
as explained in the text. Panel b: the change of the measured 
angle 7 as a function of the number of steps k away from the 
deepest particle. The angle settles on a value that depends 
only weakly on k. 

belong to the critical set. 



VI. CONCLUDING REMARKS 

In conclusion, we have used the method of iterated con- 
formal maps to compute accurately the harmonic mea- 
sure of DLA clusers of moderate size. We have explained 
that we must use the full power of the method in order to 
overcome the strong contraction of the regions on the unit 
circle that belong to the deep fjords. By iterating back 
and forth, using the fact that we own the history of the it- 
eration scheme, we could resolve probabilities as small as 
1CP 35 . Using this data we could establish beyond doubt 
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maps the outside of the unit circle to the inside of a 
wedge with opening angle 7 C = it /a, where a > 1 allows 
7 C be vary between and tt. 

To calculate the electric field -E we need the inverse of X' 



z a +i 



z™ - i 

From here we see that x _1 (0) = -1 and 



-1 ( p< Ml±l)/2a) _^ 1 



as p 



(A2) 



(A3) 



Thus the unit circle is unfolded onto the wedge; shifting 
the point w = — 1 to the origin z = 0; and rotating and 
stretching the upper half circle onto the real axis and the 
lower half circle to the other ray of the wedge p e m l a . 
The electric field follows from its definition: 



FIG. 8: The distribution of angles 7 as determined by the 
procedure exemplified in FigJ^ over the set shown in Fig. ^. 



that the generalized dimensions (or, equivlently, the f(a) 
function) lose analyticity at a negative value of q, imply- 
ing the existence of a max . In order to understand the loss 
of analyticity we offer a geometric picture. We identified 
the critical set on the cluster as having harmonic proba- 
bilities that belong to the power law tail of p(X n ). Con- 
sidering this set we identified fjords that can be modeled 
as wedges of characteristic angle. Taking such wedges 
as a model for the fjords of the critical set, we found a 
value of q* which is very close to the one computed us- 
ing other methods. We thus propose that the point of 
non-analyticity can be interpreted as resulting from the 
power-law dependence of the harmonic measure in the 
fjords belonging to the critical set. 



E(z) = 



— log(x 



2a 



1 



\z\ \z a + z~ a \ 



(A4) 



On the real axis close to the center z = 0, i.e. for 
z = p <C 1 it becomes 



E{P) = ^ 



1 



p p a + p-' 
while for large p is goes like 



2ap° 



for p < 1 (A5) 



, , 2a 
E(p) = — 



1 



2a p 



-(a+l) 



p p a +p- 

Exactly the same relations hold for the upper ray. 



for p > 1 
(A6) 



2. The Probability Measure for A„ 
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APPENDIX A: WEDGE MODEL FOR THE 
FJORDS OF THE CRITICAL SET 



The linear size of the particles in mathematical space 
/\n is proportional to the electric field. 



(A7) 



\ n (0) =V\ E($ (e i9 )) 



Thus the probability measure of the A„ is directly related 
to the probability measure of the electric field. Since E 
is the same for the two rays of the wedge, it is sufficient 
to consider it only on the real axis. Starting from the 
uniform distribution of the 0-values in mathematical 
space, it follows: 



1. The Conformal Map and the Electric Field 

The conformal function 



X{w) 



1 



w 



l/a 



(Al) 



dP (6») dP(E) 
dE 



dO 




dP 


dE 


dE 


dp 


dP 


dE 


dE 


dp 



dE 



dp 



dp 
d0 



^-log^-Hp)) 



E~ 



(A8) 
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or 



dP{E) 



dE 
dp 



E dE 



(A9) 



The derivative of the electric field follows from (A5) or 

(M) 



dP [y/Xn)~y/Xn (A13) 



dE , 
dp 

Yielding 



(p) 



2a \(a+ l)p a - (a- l)p~ a \ 
P {p a + p- a f 



(A10) 



For the probability density of A„ this means 



dP(E) 



{a+l)p a -(a-l)p- a \ y ' 



dP(Xn) = 



dP (VK) d\. 



For small p corresponding to a small field and thus small 
VAT we get: 



d-\/Xn 2 \fX~^ 

- a ~ 2 
X n d\ n 



(A14) 
(A15) 



dP(E) ~ pdE ~ E< 



(A12) 
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